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o : 

O . Abstract 
(N ■ 

Q I CP and T asymmetries in neutrino oscillations, in vacuum as well as in matter, are 

■ expressed in terms of invariant functions of lepton mass matrices. 

o : 

cN 1 Introduction 



J> ' Four decades have passed since the unexpected discovery of CP violation, originally seen in two 
^ ■ decay modes of ■ For about 25 years the superweak ansats by Wolfenstein |2j accounted for all 
psj observed effects, including CP violation in the semileptonic decay modes of Kl, and provided a 
C^l simple and intuitive explanation for why CP violation was not seen elsewhere in particle physics. 
^ [ Nowadays the standard framework for understanding CP violation is the electroweak model [Sj 
O I (hereafter referred to as the standard model) and the Kobayashi-Maskawa scheme |0] within it. 
^ • Indeed the remarkable experimental determination of the quantity e'/e in decays and the 
^ discovery of CP violation in the decays of B-mesons, at SLAC and KEK are in agreement with the 
P-i, predictions of the standard model with three families. Nonetheless, understanding CP violation 
r-| still remains a great challenge, within a broad area in physics. One faces deep questions such 
' as what is the source of CP violation that goes into generating the baryon asymmetry of the 
universe and why is the theta parameter of QCD so small that it has not been seen. There are 
also some perhap simpler questions, one being: is there CP violation in the leptonic sector? The 
minimal standard model gives a no as the answer to the latter question because it assumes that the 
neutrinos are massless. But nowadays there is evidence for neutrino oscillations, a phenomenon 
that requires massive neutrinos Therefore the minimal standard model needs to be modified, 
but we don't know how. In this article, 1 shall address the issue of CP violation in neutrino 
oscillations using the method of invariant functions of mass matrices ([H|,[7J). In this article, first 
a short introduction to this method is presented in the next section followed by application to 
neutrino oscillations, in vacuum as well as in matter, in the following sections. 
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2 Invariant functions of mass matrices 

Consider first the quark sector of the standard model with three families. The identity of the 
quarks is encoded in the three-by-three quark mass matrices M„ and Md, for the up-type and 
down-type quarks respectively. However, these mass matrices are basis dependent. Given any 
pair Mu, one can obtain other pairs through unitary rotations, as will be described below, 
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without affecting tlie physics. The measurable quantities must be basis independent and therefore 
they are "invariant functions" under such rotations. These functions were introduced in [Hj and 
studied in detail in [7]. Actually, what enters, in the standard model, is the pair 

Su = MuMl , Sd = MdMl (1) 

For simplicity, we shall refer to these quantities as mass matrices. This should cause no confusion 
because the underlying mass matrices, M„ and Md, do not enter in what follows. 

An invariant function f{Su, Sd) is a function that does not change under the transformation 

Su ^ XS^X\ Sd XSdX^ (2) 

where X is an arbitary three-by-three unitary matrix. Evidently, the traces of powers of the above 
quantities, tr[S^ S"^), are such invariant functions. Note that the corresponding determinants are 
not independent invariant functions because any determinant can be expressed as a function of 
traces. For a detailed discussion of these invariants see pj. 

When dealing with CP violation in the standard model with three families, a central role is 
played by the commutator of the quark mass matrices, [Su, Sd\ (see (0, [H])) The determinant of 
this commutator is an invariant function of mass matrices given by ([H], [H]) 

det [Su, Sd] = 2iJ v{Su)v{Sd) (3) 

where J is the CP-invariant of the quark mixing matrix V, 

ImiV^jVpkV^kV^^) = J 5^e«;3^e,fci (4) 

J is equal jH] to twice the area of any of the six by now well-known unitarity triangles. The quan- 
tities v{Su) and v{Sd) are Vandermonde determinants as follows. Denoting the three eigenvalues 
of Su by Xi, Xi = m\, X2 = m^, x^, = rri^, we have 

viSu) = ^ijkXjxl = {ml - ml){ml - ml){ml - ml) (5) 

and similarly v{Sd) = {m'^ — m^)(m^ — ml){m1 — m^). 

In the standard model with three families, the nonvanishing of the above determinant gives 
the if and only if condition for CP violation in the quark sector. In fact it manifestly unifies the 
14 conditions needed, such as the condition that no two quarks with the same charge are allowed 
to be degenerate if CP is to be violated, or the conditions that none of the mixing angles nor the 
phase angle is allowed to assume its maximal or mininmal value. In other words, the commutator 
automatically keeps track of these requirements and is thus useful for checking whether a specific 
model violates CP or not (for a pedagogical discussion see [TUj). 

The determinant in Eq. appears naturally in computations involving CP violation when all 
the six quarks enter on equal footing. Examples are the renormalization of the ^-parameter of 
QCD by the electroweak interactions and the calculation of the baryon asymmetry of the universe 



2 



in the standard model. The determinant, in a truncated form, enters in many more computations 
such as when computing the electric dipole moment of a quark, say the down quark. Since in 
such a computation, the down quark appears in the external legs, it is tacitly assumed that we 
know the identity of this quark, i.e., 7^ and 7^ rrih. Therefore the factors (m^ — m^) and 
{ml — ml) will be missing but all the other factors will be present. 

The above commutator is the simplest in a family of comutators of functions of mass matrices 

C{f,g)^[fiS^),g{S,)] (6) 

/ and g being functions that are diagonalised with the same unitary matrices that diagonalise 5*^ 
and Sd respectively. The determinants of these commutators, which are also invariant functions, 
are given by 

det[f,g]=2zJv{f)v{g) (7) 
where f (/) = v{f{Su)), and v{g) = v{g{Sd)). More explicitly 

V{f) =Y.'^,^f^fl = - /2)(/2 - /3)(/3 - /l) (8) 

The fj denote the three eigenvalues of the matrix f{Su) and the quantities related to the down- type 
quarks are defined similarly. 

An essential point is that Eq.(|7)) holds irrespectively of whether / and g are hermitian or not. 
This property makes the above formalism applicable to neutrino oscillations, as we shall see here 
below. 



3 Neutrino oscillations in vacuum 

With massive neutrinos, CP violation in neutrino oscillations could manifest itself by the difference 
of the rates of reactions z/^ — > z/^ and — > up. Here a and /3 (a 7^ /5) stand for e,fi,T. The 
difference of these rates is found to be proportional to 

sm(20i2) + sm(2023) + sm(203i) = -4sm(0i2)sm(023)sm(03i) (9) 

0,, = (m| - m2)e, e = ^ (10) 

m's being the neutrino masses; E denotes the neutrino energy and L is the distance from the source 
(for a recent review see |I2])- Because of CPT invariance, the above combination of sm's also 
appears when testing time-reversal asymmetry by comparing the rates of — I'p and vp ^ v^. 

As was noted before [7j, in models with quark- lepton universality the results for the quark 
sector, presented in the previous section, can be extended to the leptonic sector by trivial sub- 
stitutions, M„ — i> My and Mi where M^ and Mi denote the neutrino and the charged 
lepton mass matrices. An interesting question is then: how is the sum (or product) of the above 
three sin's related to the commutator of lepton mass matrices? One would expect that there be a 
relationship because, as far as neutrino oscillations are concerned, the leptonic sector is essentially 
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a copy of the quark sector. An essential point is that the possible Majorana nature of neutrinos is 
known not to be important in the computation of the rates (see ^M)- The relevant formulae 
are obtained by assuming that there are three active neutrinos and that their "effective" mass 
matrix is three by three. This pattern emerges in many models, some based on the see-saw mech- 
anism (for a recent review see ^Hj ) as well as in models a la Weinberg [IGj where the effective 
operator that generates the neutrino masses involves only the three active left-handed neutrinos. 

To answer the question concerning the relationship between the above sum or product of the 
three sin's and the lepton mass matrices, we introduce, in analogy with the case of the quarks, 

= M^Ml , Si = MiM} (11) 

where My and Mi are the three-by-three neutrino and charged lepton mass matrices respectively. 
We introduce the commutators 

A± = [e^^'^'^ , Si] (12) 

where the unitary matrices = e^^*^"^" are inverses of one another and A"*" = (A")"!". The 
determinant of these commutators are invariant functions of lepton mass matrices. Using Eq.Q 
we have 

detA^ = 2i Jy v{Si)v{e^^'^''^'') (13) 

Here Ji, is the leptonic analogue of the CP invariant of the quark mixing matrix. Here it is more 
convenient to use the index u instead of "lep" (for leptons) because when dealing with oscillations 
in matter, in the next section, the notation is easily generalised. and Jp will then denote the 
corresponding quantities for neutrino and antineutrino oscillations in matter. 

Furthermore, just as in the case of the quarks, Ji, is simply twice the area of any of the six 
leptonic unitarity triangles. The two Vandermonde determinants in the above equation are given 
by 

v{Si) = {ml - ml){ml - ml){ml - ml) (14) 

and 

t,(e*2iC5.^) ^ ^8ie^^^^*'^^^sm(0i2)sm(023)sm(03i) (15) 

where trS^, = {mf + m"^ + m"^) and the m's are again the neutrino masses. The exponential factors 
are the determinants of U^. Putting these results together, we find 

detA^ = ±16v{Si) [J^.sm(0i2)sm(023)sm(03i)] e^^^^*"^" (16) 

Note that the phase of the determinant is determined by the sum of the neutrino masses. Evidently 
any of the relations detA^ ^ provides the necessary and sufficient condition for having CP 
violation in reactions ^q, — >■ and Ua ^p- The presence of the factors involving the masses 
of charged leptons is essential. These keep track of the identity of neutrinos. For example, for 
me = the electron and muon neutrinos would be indistinguishable and therefore there would 
be no CP violation. detA'^ ^ also provides the if and only if condition for T-violation when 
comparing the reactions Va and Va- 
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4 Neutrino oscillations in matter 



At a first glance, the formalism for CP and time-reversal violation in neutrino oscillations in matter 
looks deceptively similar to that in vacuum. It would seem that all one needs to do is to distinguish 
the masses and mixings in matter by simply putting say primes on corresponding quantities in 
vacuum, as we shall do here below, and by introducing explicit indices v and v to keep track of 
whether we are dealing with neutrino or antineutrino propagation in matter. However, there are 
subtle points to be taken into account, as we shall see soon. 

Consider first the case of neutrinos. To the leading order, in a frame where the charged lepton 
mass matrix is diagonal, the neutrino mass matrix, Sy = M^Ml, is replaced by an effective mass 
matrix in matter, S'^ = M'^M'J , given bv ([T7j. |18j) 

Sl = Sy + n, n = diag{p, 0, 0) (17) 

where p is proportional to the neutrino momentum as well as the density of the electrons in matter 
(for a review see, for example, ^H]). Here, we shall treat p as a constant. The reality is generally 
far more complicated. But here we are primarily interested in exploring the structure of matter 
oscillations rather than making realistic calculations. Thus we construct the matter analogue of 
any of the two vacuum commutators and take its determinant. We take A+, drop the superscript 
and define 

^ p^f^^ (18) 

From Eq. ()13p we have 

detAl = 16 [Jl v{Si)sin{(f)[^)stn{(j)'^^)sin{(f)'^^)] e^^^*"^^ (19) 

Here the primed quantities are the matter analogues of the unprimed ones in the previous section 
and, as before, the phase of this determinant is fixed by the sum of the neutrino masses, now taken 
in matter. Thus trSl = (mf + m2 + rnf); m"s being the neutrino masses in matter, etc. The 
subscript u is a. reminder that we are dealing with neutrinos propagating in matter. Evidently 
Jl ^ Jy as p ^ 0. 

Actually, this determinant is not what enters when testing CP violation in neutrino oscillations 
in matter through the rates discussed above. It will, however, enter if we were to do the science- 
fictional experiment of comparing the rates of z/q, —>■ vp in matter with that of Va ~^ in the 
corresponding antimatter. It is ironic that the very same baryon asymmetry of the universe, 
that owes its existence to CP violation, forbids us to test CP violation in matter in the above 
"straight-forward" fashion. 

Turning now to antineutrino oscillations in matter, the effective mass matrix, = M^Mp , to 
the leading order is given by 

S'p = Sy-U = Sy- diag{p, 0, 0) (20) 

Because of the different sign of the added term for neutrinos and antineutrinos, the effective 
neutrino mass of the j-th neutrino, m'j, is in general not equal to that of the corresponding 
antineutrino, m' . Nor are the mixing angles and thus the corresponding J's in general the same. 
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This is the reason why we need to introduce not only primes to indicate the presence of matter 
but also appropriate v and v indices to keep track of whether we are talking about neutrinos or 
antineutrinos. Defining the commutator relevant to antineutrinos, we have 

rfetA; = \^v{Si) {j'^s%n{^^^)sinQ^^.^s%n{^^^\ e^^^*"^^ (21) 

Depending on the density profile of matter, the above expressions could, however, be relevant 
when testing time-reversal violation. The quantity detN^ would enter when comparing the rates 
of Va — ^ 1^/3 and vp Va in matter and detA'^ for corresponding antineutrinos processes again in 
matter. 

Since the matter contribution to the mass matrices, 11, commutes with the lepton mass matrix 
Si (which is diagonal), we have 

[S,,Si] = [Sl,Si] = [S',,Si\ (22) 

Taking the determinants of these matrices and removing the common factor involving the charged 
leptons, gives 

J. v{S,) = Jl viSl) = 4 viS',) (23) 

This result can be found in Refs.|,2.Qj and [21.^. Again v{S^) = {mf — m\){m\ — m|)(m| — m^), 
the m^'s being the neutrino masses in vacuum; v{Sl) and f (S'p) are the corresponding expressions 
with neutrino and antineutrino masses in matter respectively. A more detailed discussion of some 
further matter effects is given in Refs. 0| and 

We shall now consider the next simplest commutator in the series, i.e., [Sj, Si\. The deter- 
minant of this commutator can be written down immediately, in terms of the matter quantities, 
using the general formula Eq.((Z|). We have 

det[S';},Si] =2iJlv{Si)v{S';}) (24) 

where 

v{Sl') = wiS'MSl) (25) 

wiSl) = (mf + m'i)im'i + mf)imf + mf ) = lutrSlf - trS'^} (26) 

The determinant in Eq.(j^J= can also be computed directly, in terms of vacuum quantities, by 
substituting SI = + H- A straight-forward calculation gives 

det [Sl^ Si] = 2iJ, v{Si)v{S,) [w{S,) + ph + p'h] (27) 

where 

h = {trS,f-Y,\yc^i?< 

a 

h = trSu - ^ I Vai\^m\ 

a 
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Here all the masses and mixings refer to vacuum quantities. Note that the coefficient of vanishes 
because the 11 [S*", 5*/] 11 = 0. Equating the two expressions for the determinant and using that 
J'u ^i^l) = Jy v{Sy) gives 

w{Sl)=w{S,)+ph+p^h 
Rewritting this relation in terms of traces gives 

[trSlf - trS'^ = {trSyf - trS^ + ?>p{trSyf - 3tr(5^n) + ?,phrSy - ?,tr{SyIl^) 

This relation can easily be checked by computing its left-hand side, in terms of vacuum quantities, 
using again S*^ = 5*1^ + 11. 

For antineutrinos the corresponding results are easily obtained by flipping the sign of 11 and 
therefore also that of p. 



5 Oscillations in matter with low density 

The results obtained in the previous section hold to all orders in the matter-related parameter p. 
Here we would like to examine how the results look like in the low density limit, i.e., when terms 
of order p^ and higher can be neglected. We shall consider the case of neutrinos, extension to 
antineutrinos being trivial. Beginning with the matter commutator in Eq. ()18|l we write 

detM^ = detAy + pR+ 0{p^) (28) 

where the flrst term in the RHS is the vacuum contribution. To compute R we expand the 
exponential in Eq. (fTH|l using Eq. (|T7jl and flnd 

= + E ^ E [^r'^'ns^ s,] + o{p') (29) 

n=l ■ fc=0 

Taking the determinant of the RHS we flnd 



R=-v{Si) E{a,P,j,a)F{a,P,j,a) (30) 

where 

E(a,/3,7,a) = J] e,,,(al)^(«r)(/3r)*(/?s)(7s)^(7t)(at)*(al) (31) 

rst 

F{a, (3, 7, ^) = ( , % )e''am,+^.) (32) 

Here, for simplicity, we have introduced the short-hand notation (ar) = Var, V being the vacuum 
lepton mixing matrix. Note that all the dummy indices in the above sums run from one to three. 
Note that E*{a, (3, 7, cr) = —E{a, 7, f5, a) and since F is symmetric under 7 as well as under 
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a ^ o" we only need the corresponding symmetric part of the function E. This quantity is found 
to be 

E{a,P,'J,a) -2iJ^er,ap{5aaSl3y + (1 - |^7l|^)(l - Saa) " (1 " |Kl 1^)^/37} (33) 

where the dummy index rj is summed from one to three. Multiplying the above factors, in Eq. fj3Uj) . 
yields 

pR = -WiCr]v{Si) J^Xe^'^''^" (34) 



where 



X = - {(l)af3COs{(f)al3) + 2 1 | ^ [cos(0^^) Sm(0^a) - (f)al3COs{(f)af3)]} (35) 



Note that X vanishes, as it should, if any two neutrinos would be degenerate. Furthermore, the 
measurable quantities \Vai\ are themselves invariant functions of lepton mass matrices. They can 
be extracted from the mass matrices with help of projection operators P. In our case we have 



.2 lwD.c^TT^_ ltr[{S,-ml){S,-m'^)U] 



Kir = -tr{P^{S,)U) = ^ (36) 

where, as before, ^(m^ — m^) = 0q,^. The identity 

Sin{(j)ap)c0s{(j)prj)sin{(j)r,a) = ^[cOs(20„/3) + COs(20^^) + C0s{2(j)r,a) - 1] (37) 

vahd for (j)a/3 + (I>f3r} + 4>r]a = 0) cau be used to rewrite X noting that the RHS of Eq. ljHTj) has the 
same value for every term in X. 

Returning to Eq. (j28|) . the LHS is a function of matter variables and the RHS contains only 
vacuum variables. Furthermore, the vacuum and the first order correction terms have the same 
phase. Removing common factors, we find 

Jl sm(0i2)sm(023)sm(03i)e^*^^ = Jiy{sm(0i2)sm(023)sm(03i) - i^pX + 0(p^)} (38) 

Thus 

Jl sin{(j)[2)sin{(f)'2s)sin{(j)'^^) = Jy [sm(0i2)sm(023)sm(03i) + 0{p^) (39) 
Substituting the relation J'^ v{Sl) = v{S,,), from Eq. (jHI?|) . we find 

— ^sm(0;2)sm(023)sm(03J = ^^sm(0i2)sm(023)sm(03i) + (9(p^) (40) 
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